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DISCOVERY OF ry-MESIC NUCLEI 

G.A. SokoP, L.N. Pavlyuchenko 
P.N. Lebedev Physical Institute, Moscow, Russia 

Abstract. First experimental results from a photoreaction which confirm an ex- 
istence of eta-mesic nuclei are discussed. The experiment was performed at 
the 1 GeV electron synchrotron of LPI. Two values of the end-point energy of 
the bremsstrahlung beam were used, 650 and 850 MeV (below and above rj- 
photoproduction threshold on the nucleon). Correlated 7r+n pairs flying from 
the ^■^C target were detected by two TOF scintillator spectrometers placed in op- 
posite directions and transversely to the 7 beam. Measured energy spectra provide 
a strong evidence for a creation of eta-mesic nuclei as intermediate states of the 
reaction studied. 

1 Introduction 

A modern view of the structure of atomic nuclei as a system of protons and 
neutrons appeared in 30s years of the 20th century after the discovery of the 
neutron. Since then many new particles and their resonance states were found 
which led in 60s to a development of the conception of quarks. Accordingly, 
a question arose of a possible existence of nuclear bound systems including, 
apart from protons and neutrons (nucleons), new particles too. 

In 50s years, hypernuclei — a new sort of nuclei consisted of A or S hyperons 
(apart from the nucleons) — were discovered. The A and S hyperons are close 
analogues of nucleons as for their masses and the quark content. Their masses 
are only bigger by ~ 200 MeV than the nucleon mass (~ 940 MeV) and they 
consist of 3 constituent quarks as the nucleons do, one of them being however 
the strange (s) quark. 

2 rjN interaction and 77-nuclei 

An idea that a bound state of the 77-meson and a nucleus (the so-called 77- 
mesic nucleus) can exist in Nature was put forward long ago by Peng [1] . The 
rj-meson does not have an open strangeness. It consists of 2 quarks (more 
exactly, a quark and antiquark with the total isospin 0) of different flavors, of 
which Ri 50% is an ss pair. The mass of 77 is 547.5 MeV, i.e. about 1/2 of the 
nucleon mass. The suggestion of J.C. Peng was based on the first estimate of 
the rjN scattering length a^Ar, 

a^AT = (0.27 -f i • 0.22) fm, (1) 

derived by Bhalerao and Liu [2] from a coupled-channel analysis of the reac- 
tions vrA^ ttN, ttN rjN and ttN ttttN. Owing to Rea,,Ar > 0, there 
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is an average attractive s-wave potential between a slow rj and a nuclcon. For 
extended nuclei, such an attraction should be sufficient for making the j^-meson 
bound, provided the life-time of the 77-meson in the nucleus is not too short. 
A quantum-mechanical consideration done by Liu and Haider [3] and based on 
the rjN potential corresponding to Eq. (1) predicted that bound states of the 
ry-meson and a nucleus A exist for all j4 > 11. Later on, this conclusion was 
strengthen. More sophisticated coupled-channel analysis [4,5] taking into ac- 
count both resonance and nonresonance contributions arrived to the scattering 
length Re a,, AT which is about 3 times bigger than the very first estimate (1): 

a^jv = (0.75 + i • 0.29) fm [5]. (2) 

For such a large arjN, rj-mesic nuclei are predicted to exist for all nuclei with 
A > A. With slightly larger o^at, r/-bound states are to be possible for A = 3 
and even for A = 2 [6] . It is worth emphasizing that elementary rjN scattering 
amplitudes are theoretically derived from other reactions like ttN — > TriV and 
ttN — > r/N through extrapolations based on a factorization [2] or, in the latest 
works, on unitarity constraints [4,5]. Since, however, not all important chan- 
nels are involved into those extrapolations (missing are, e.g., t:N KA and 
ijN KA), it is not clear how reliable are the obtained results. A difference 
between Eqs. (1) and (2) may give a hint about real uncertainties. There- 
fore, experimental studies of bound states of various rjA systems would greatly 
contribute to learning elementary t]N scattering. 
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Figure 1: Energy dependence of the amplitude /^jv for the process rjN — > rjN [5]. 

Real part of the rjN scattering amplitude frjN , the threshold value of which 
is just equal to the scattering length a,,Ar, remains positive up to the kinetic 
energies of rj below 70 MeV [5] (Fig. 1). This means that an effective rjA at- 
traction exists in a wide near-threshold energy region E^™ « 70 MeV. The 
attractive forces in the final state have to lead to a near-threshold enhance- 
ment in the total and differential cross section of real-77 production by different 
beams. Such an enhancement was indeed observed in several reactions includ- 
ing p((i, ^He)?] [7,8] and d{d,^lle)ri [9,10], thus supporting the existence of the 
rjA attraction even for the lightest nuclei. Nevertheless, all these experiments 
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which have deal with rj in the final state cannot directly prove that bound r]A 

states do really exist. A well-known counter example is provided by the NN 
system in the ^Sq state which has a virtual rather than a real level and it has 
a negative rather than positive scattering length. 

3 A previous experience 

Two attempts to find Tj-nuclei in the missing mass spectrum of the reaction 
7T+A pX were performed at BNL [11] and LAMPF [12] soon after the first 
theoretical suggestions [1,2]. Both the experiments failed to find a signal of the 
?7-nuclei, perhaps owing to their much bigger total widths than then expected. 

A new experiment started in 1994 at the 1 GeV electron synchrotron of 
Lebedev Physical Institute aiming at a search for ry-nucleus bound states in 
photoreactions. In 1998, the first experimental data were obtained [13] which 
indicated that a bound state of the 77-meson and a nucleus ^ = 11 does exist. 

4 LPI experiment 

Below we describe in some details a method used in [13] for identification 
of the bound states of the r?- meson and a nucleus, as well as a procedure of 
measurements and data analysis. 

As a preliminary remark note that the use of the photon beam for formating 
//-nuclei has certain advantages when compared with the use of the pion beam. 
The reason is that 

- photons interact with nucleons throughout the whole nuclear volume whereas 
pions strongly interact only with nucleons near the nuclear surface; 

- high photon flux of order 10^ — 10^° s^^ in the energy range of AE^ ~ 200 
MeV provides quite an adequate reaction rate. 

5 Method of identification of 77-mesic nuclei 

An important idea of the work [13] was in searching for decay products of 77- 
nuclei which are ttN pairs arising from decays of the S'ii(1535) resonance inside 
a nucleus [14]. A mechanism of r;-nuclei production and decay in photoreac- 
tions is schematically shown on Fig. 2. At the first stage, the incoming photon 
produces a slow f/-nic!son and a fast nucleon, the latter leaving the nucleus. 
Then, in the nuclear medium, the ?7-meson undergoes multiple elastic scatter- 
ing, r]N — > 5ii r}N — > 5ii At the last stage, the jy-meson annihilates 

after interaction with a nucleon, 77 TV ^ 5*11 ^ tiN . so that a tiN pair emerges 
which leaves the nucleus too. The 511(1535) resonance plays a fundamental 
role in this dynamics. It provides production and annihilation of the r/-meson 
and it also predetermines an attractive ry-nucleus interaction. The chain of 
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Figure 2: Mechanism of creation of a ?j-mesic nucleus and its decay into a -kN pair in a 

photoreaction. 

elastic processes r]N — > Su -qN . . . Sn in the nuclear medium results in 
randomizating momenta of the 5*11(1535) resonance and the r;- meson. Fermi 
motion of nucleons is important for this randomization because the 77-meson 
scatters on nucleons having random velocities and directions. The most im- 
portant result of the randomization is the isotropic momentum distribution of 
decaying 5*11(1535) resonances. It leads to an isotropic distribution of the ttN 
pairs and, since the momentum of 5*ii(1535) in the nucleus is small, to the 
opening angle Oj^n close to 180°. This is a basic idea of the method of ?7-nuclei 
identification suggested in [14]. The method consists in detection and energy 
measurements of components of wN pairs from 5*ii -resonance decays in the 
nuclear medium. It is worth mentioning that Tr+n pairs flying transversely to 
the photon beam cannot be produced via the one-step reaction 7p — » 7r+n in 
the nucleus when the photon energy is as high as 700—800 MeV. 

6 Simulation of ttN events 

Theoretical estimates [13,15] show that the binding effects lead to a full dom- 
inance of the reaction mechanism related with multiple rcscattcring of 77 and 
with a formation of the intermediate ?7-nucleus over a mechanism of non- 
resonance (background) production of the wN pairs in the subthreshold region 
of the invariant mass y^i^Jwv < 'tijj + itin for the subprocess rj + N ^ tt + N . 
A peak in the mass distribution of ttN is theoretically expected in this region. 
This is illustrated in Fig. 3 where the spectral function S{E) of the (kinetic) 
energy i? of 77 in the nucleus is shown which is proportional to the number of 
rjN collisions that the 77 undergoes when travels through the nucleus. Another 
spectral function, S{E, q), shows a distribution of the produced ttA'' pairs over 
their total energy E + + niN and the total momentum q. A presence of 
the rjN attractive potential U produces strong enhancements in these spec- 
tral functions in the energy-momentum region corresponding to the bound rjA 
states. 
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Figure 3: Spectral functions S{E) and S{E, q) (in arbitrary units) of tlie kinetic energy 
E and momentum g of in the nucleus. They are found with a rectangular-well optical 
potential simulating the nucleus ^^C. For a comparison results obtained with dropping out 
the attractive (i.e. real) part of the rjA potential U are also shown. 




Figure 4: Layout of the experimental setup. Shown also are the time-of-flight spectra in the 
TT (left) and n (right) spectrometers. 

7 Experimental setup and procedure of measurements 

The reaction studied in our work was 

7 + ^ + iiB (iiC) ^7r++n + X, (3) 

in which energies and momenta of of tt"*" and n are correlated. Among four 
different TriV channels of S'ii(1535) decays we have chosen the Tr+n channel 
because: 

- detection of the Tr^-decay modes is more expensive since this requires a use 
of the 47r geometry of the 7-spectrometer; 

- detection of the tt~p pairs with proton energies of about 100 MeV dictates 
a use of rather thin targets and leads therefore to a reduced count rate. 

At the same time the selected decay channel (Tr+n) allows the use of thick 
targets and a rather simple setup. Average energies of nn pair's components are 
(E^) = 300 MeV and (£'„) = 100 MeV. The setup used included two arms of 
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time-of-flight (TOF) scintillator spectrometers (Fig. 4). The pion spectrometer 

contained two scintillator blocks of 50 x 50 x 2 cm"^ and 50 x 50 x 5 cm^. The 
neutron spectrometer contained an anti-coincidence scintillator detector of the 
size 50 X 50 X 2 cm^ and four scintillator detectors of the size 50 x 50 x 10 
cm"^. Efficiency of the anti-coincidence (veto) detector was 90%. Efficiency 
of the neutron detecting was about 8 10%. Each scintillator detector had 
four photo-tubes FEU-63 placed on the edges of the scintillator block; they 
were used for determination of the coordinates of particles passed through the 
detector using time differences of light signals coming to the photo-tubes. An 
accuracy of the coordinate determination was about Aa; = A?/ = 1.5 cm. An 
accuracy of the determination of time was At = 250 ps. Measurements in 




Figure 5: Raw Tr+n-event distributions over the particle velocities /3„ and f3n for the "cali- 
bration" (a), "background" (fe) and "effect -|-background" (c) runs. 



the experiment were made in three different runs. In the "calibration" run 
(a) both 7r+ and n spectrometers were positioned at angles (6*^) = +50° and 
i^n) = —50° with respect to the photon beam direction and the bremsstrahlung 
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Figure 6: Unfolded twodimentional distributions over the velocities /3's of the 7r+n events 
with the end-point energy of the bremsstrahlung spectrum fi-^max = 850 and 650 MeV. 

beam end-energy E'^max = 650 MeV was used. This configuration is suitable for 

measuring quasi-free singlc-pion photoproduction jp — > Tr~^n on protons in the 
nuclear target. It was used for final tests of both arms of the setup since energies 
of the detected particles were close to expected energies of 7r+n pairs arising 
from the 5*11(1535) decays. In the "background" run (&) the spectrometers 
were shifted to (6*^) = +90° and (6'„) = —90° and the beam energy was still 
Ejma^ = 650 MeV. In such a configuration single-pion production events were 
eliminated and most of the detected events were related to a quasi-free doublc- 
pion photoproduction ■^N ttttN. In the last, "effect -l-background" run (c) 
the configuration of the setup was as in the "background" run (6) with the 
except for the end-point energy which was increased to E'-ymax = 850 MeV thus 
making rj production possible. 

8 Data analysis and results 

Fig. 5 shows a two-dimensional spectrum of the detected particle's velocities /3 
for all three runs. The measured velocities 0i = Li/cU are subject to fluctu- 
ations stemming from errors Sti and SLi in the timc-of-flight t, and the flight 
base Li measured in the experiment. Such fluctuations are clearly seen in the 
case when tt'^tt'^ (actually photons or relativistic electrons or positrons) hit scin- 
tillators. Therefore, experimentally observcxl velocities are close but not equal 
to 1 (see in Fig. 5). Therefore, an experimental /3-resolution of the setup can 
be directly inferred from the ir^n^ evens. Then, using this information and 
applying a statistical method of solving the inverse problem described in [16], 
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Figure 7: Distribution over the total ki- 
netic energy of the 7r+n pairs for the 
"effect+background" run (left panel) and for 
the "background" run (right panel) obtained 
after unfolding raw spectra. 



Figure 8: Distribution over the total 
energy of the tt+ti pairs after a sub- 
traction of the background. A distri- 
bution corresponding to a decay of a 
free 5ii(1535) resonance is shown for 
a comparison. 



one can unfold the experimental spectrum, obtain a smooth velocity distribu- 
tion in the physical region Pi < 1 (Fig. 6) and eventually find a distribution of 
the particle's kinetic energies Ei — Mi[{l — Pf)^^^^ — 1]. Finding Ei, we intro- 
duced corrections related with average energy losses of particles in absorbers 
and in the detector matter. Fig. 7 shows a two-dimensional energy distribution 
of TT^n pairs obtained through the above data evaluation for two runs at the 
beam energies of 650 and 850 MeV. At 650 MeV, Figure 7 shows a smooth 
decrease of the number of 7r+n pairs with the total energy of the pair. At 
850 MeV, a resonance structure is observed in the kinematical region of the 
511(1535) resonance. We believe that such a structure is directly related to 
77-nuclei formation and decay. 

Of the most interest is the distribution of the 7r+n events over their total 
energy E'tot = En + E^^. Subtracting a smooth background (Fig. 7), we have 
found a 1-dimensional energy distribution of the 7r+ri events presumably coming 
from bound r;- mesons decaying in the nucleus (Fig. 8). The experimental width 
of this distribution is about 150 MeV. Its center lies by Ai; « 90 ± 15 MeV 
below the position of the 511(1535) resonance and even below the threshold 
energy m,, + ttiat = 1485 MeV, thus indicating a presence of binding effects, 
i.e. a formation of ?7-mesic nuclei. 

9 Prospects 

Studies of ?7-mesic nuclei which lie at the intersection of the nuclear physics and 
physics of hadrons promise us to bring a new information important for both 
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the fields. 77-mesic nuclei provide a unique possibility to learn interactions of 
Ty-mesons with nucleons and nucleon resonances in the nuclear matter. Data on 
the self-energy of S'ii(1535) in the medium, being interpreted in the framework 
of the chiral-symmetry models, can shed more light on the problem of masses 
of free and bound hadrons [17]. 
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